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ABSTRACT 

In this Letter, we present mean optical+NIR color gradient estimates for 5080 early-type galaxies 
(ETGs) in the grizYJHK wavebands of the Sloan Digital Sky Survey (SDSS) plus UKIRT Infrared 
Deep Sky Survey (UKIDSS). The color gradient is estimated as the logarithmic slope of the radial 
color profile in ETGs. With such a large sample size, we study the variation of the mean color gradient 
as a function of waveband with unprecedented accuracy. We find that (i) color gradients are mainly 
due, on average, to a metallicity variation of about — 0.4dex per decade in galaxy radius; and (ii) a 
small, but significant, positive age gradient is present, on average, in ETGs, with the inner stellar 
population being slightly younger, by ~ O.ldex per radial decade, than the outer one. Also, we show 
that the presence of a positive mean age gradient in ETGs, as found in the present study, implies 
their effective radius to be smaller at high z, consistent with observations. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxi es: fundamental parameters 



1. INTRODUCTION 

Since the early 1970's, observations have shown that 
the color of ETGs reddens from the outskirts to the 
galaxy center (Faber 1972). Understanding the origin of 
such color gradient can strongly cons train the scenario 
of ga laxy formation and evolution (see iLa Barbera et ahl 
l2004f l. Most studies to date have been plagued by 
the large uncertainties due to small sample sizes and 
short wavelength baselines where color gradients are de- 
rived. Despite these limitations they suggest that color 
gradients in ETGs originate from metallicity variations 
(Peletier et al. 1990; Tamura et al. 2000; Tamura & Ohta 
2003). Small, positive age gradients are also consistent 
with observations (e.g. Saglia et al. 2000.) , though their 
presence has never been revealed in the family of ETGs 
as a whole. Positive age gradients are a robust predic- 
tion of the hierarchical paradigm of galaxy formation, as 
a consequence of the centrally peaked star formation of 
dissipative merging events. 

Recent studies have shown a significant, intrinsic evolu- 
tion of ETC sizes, with high redshift galaxies being more 
compact (by ^ 50% at z ^ 1) than those of analogous 
stellar mass at z ~ 0. The origin of this variation in size is 
still a matter of deb ate and may b e the result of different 
mechanisms (see e.g. lTruiilloll2009 ') . specifically (i) the ef- 
fect of minor dry merge rs; or (ii) a puffing- up mechanism 
due to AGN feedback ()Fan et al.ll200l : or (iii) the in- 
creased amount of dissipation involved in gas-rich merg- 
ers fo rming ETGs at high-redshift (Khochfar and Silk 
[200^ . 

In this Letter, we use data from the SDSS and UKIDSS 
surveys to address the origin of the mean color gradient 
in ETGs in the nearby Universe (z<0.1), connecting it 
to the observed compactness of their high-redshift coun- 
terparts ^ . The combination of sample size, homogeneity 
of data, and large wavelength baseline, from optical to 
near- infrared (g through K bands), results in an unprece- 



dented accuracy in detecting both metallicity and age 
variations inside galaxies. In fact, while optical (SDSS) 
data are very sensitive to the effects of both metallicity 
(through hue blanketing) and age, the NIR (UKIDSS) 
wavebands are dominated by the old, quiescent stellar 
populations, allowing the effects o f age and metallicit y 
to be effectively separated (see e.g. iPeletier et al.lll99"o() . 

2. THE SAMPLE 

The sample of ETGs was defined from SDSS-DR6 with 
photometry and spectroscopy available. We selected all 
galaxies in the redshift range of 0.05 to 0.095, with r-band 
Petrosian magnitude ^ Mr< — 20. The choice of the lower 
redshift minimizes the aperture bias ( Gomez et al. 2003) , 
while the upper redshift limit makes the sample approx- 
imately volume-complete, since the value of Mr= — 20 
roughly corresponds to the apparent magnitude limit of 
the SDSS spectrosc opy (r - 17.8) at z = 0.1. Follow- 
ing iBernardi et al.l (2003), ETGs are those galaxies with 
SDSS parameters eclass < and fracDeVr > 0.8. Also, 
we selected only those galaxies with spectroscopic warn- 
ing flags set to zero, and with available velocity disper- 
sion, in the range of 70 to 420 km/s. These require- 
ments yield a sample of 39,993 ETGs. We matched 
this sample to the f ourth data release of t he UKIDSS- 
Large Area Survey (|Lawrence et al.l[2007| ). which pro- 
vides NIR photometry in the YJHK bands over a sky 
region significantly overlapping the SDSS. The matching 
was done considering only frames with the better quality 
flag {ppErrBits < 16) in all bands. For each ETG we 
selected the nearest UKIDSS detected galaxy within a 
matching radius ^ of 1". 5080 ETGs have photometry 
available in aU eight SDSS+UKIDSS fflters. 



^ We assume a ACDM cosmology with Qn 
Ho = 75 km s-i Mpc-^. 



: 0.3, Ha = 0.7, and 



^ k-corrected with kcorrectvi^l^i ijBlanton et a l. 2003'^ through 
rest-fra me filters blue-shif ted by a factor (l-|-2:o), adopting zq = 0.1 
(see e.g. IHoeg et al.ll200l) . 

The matching was done with the CrossID form orthe WFCA M 
Science Archive (see http://surveys.roe.ac.uk/wsa/index.html for 
details). Changing the matching radius to 0.5" leads to decrease 
the sample size by only five galaxies, confirming the accuracy of 
the matching procedure. 
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TABLE 1 

Statistics of Color gradients. 



Vg-X 


M 




g-r 


-0.071 ± 0.003 


0.121± 0.004 


g-i 


-0.084 ± 0.004 


0137 ± 0.003 


g-z 


-0.089 ± 0.006 


0.173 ± 0.003 


g-Y 


-0.231 ± 0.006 


0.171 ± 0.004 


g-J 


-0.279 ± 0.008 


0.204 ± 0.005 


g-H 


-0.284 ± 0.008 


0.194 ± 0.004 


g-K 


-0.297 ± 0.008 


0.219 ± 0.005 



3. COLOR GRADIENT ESTIMATES 

The grizYJHK images of each galaxy were retrieved 
from the SDSS and UKIDSS arch ives, and processed with 
2DPH0T fsee lLa Barbera et aLl l2008: hereafter LdC08). 
Galaxy images were fitted with PSF convolved Sersic 
models, allowing us to estimate structural parameters, 
i.e. the effective radius, r^, the mean surface brightness 
within that radius, < /i>e, and the Sersic index, n, ho- 
mogeneously in all the bands. The PSF was accurately 
modeled ^ in all wavebands by fitting the five closest stars 
to each galaxy with a sum of three Moffat functions, as 
detailed in LdC08. Having eight wavebands, seven dif- 
ferent color gradients can be estimated. We used the pa- 
rameters Tg, ^e, and n, to estimate color indices, g ~ X 
with X = rizYJHK, as a function of the distance, p, 
from the galaxy center. For each waveband, we estimate 
the mean surface brightness of the de-convolved Sersic 
model ^ on a set of concentric ellipses ^ , whose ellipticity 
and position angle are fixed to the value of the r-band 
Sersic fitting. The color index g — X at a given radius 
p is obtained by subtracting the mean surface brightness 
values at the corresponding ellipse. Each color profile is 
fitted in t he radial range of Pm in — OAr^ to pmax = 
(see e.g. iPeletier et al.lll990l) , by using an orthogonal 
least squares fitting procedure. We estimate the g — X 



color gradients, 
g — X profiles 



as the logarithmic slopes of the 
c(g-X)/d{\ogp). 



4. CONSTRAINING AGE AND METALLICITY VARIATIONS 

Considering the effect of age and metallicity, we write 
the equations: 



dig - X) 

dlogt 



djg-x) 
d\ogZ 



(1) 



where and Vt are the logarithmic radial gradients of 
age, t, and metallicity, Z, across the galaxy. The quanti- 
ties Xt = d{g - X)/d\ogt and Xz = d{g - X)d\ogZ ave 
the partial logarithmic derivatives g — X with respect 
to age and metallicity. The above linear equations hold 
if the color indices are continuous functions of t and Z , 
and the absolute values of V z and Vt are small. 
From the operational viewpoint, we consider two stellar 
populations and describe the properties of the galaxy at 

For the PSF fitting in all wavebands, the mean value of the re- 
duced distribution amounts to ~ 1.04, while the corresponding 
90th percentile is 1.2. 

^ Using the de-convolved models makes the procedure insensi- 
tive to seeing effects. We also verified that color gradients do not 
correlate with PSF fitting values as well as with the error on 
the seeing FWHM, as estimated by 2DPHOT from the width of 
the sure star locus (see LdC08). 

^ Ellipses are equally spaced in equivalent radius p by 0.01re,r , 
where r is the the r-band r^. 



X 

>- 

N 



o - 



IN 



II o 

^ I 

X 

I 

en 



d 

I 



"T 1 1 1 1 1 1 1 1 1 1 1 r- 



i ' ' ' ' 1 ' ' 


' ' 1 ' ' ' ' L 






" M05 , r- 


BC03 


T-IGyr 






' +^ 


I'll 1 1 


CB08 ■ 
■ . 1 



-0.45 -0.4 -0.35 -0.3 



_l I I L_ 



10"^ 1.5x10'^ 2x10^^ 



A (A) 



Fig. 1.— 

tion of the effective wavelength of the filter X. The error bars mark 
2(T standard errors. The grey region shows the set of age and metal- 
licity models that best match the observations (see text). The dots 
in the inset plots the best-fitted Vz and Vt values obtained from 
the CB03 SSP models by shifting mean color gradients according 
to their uncertainties (see text). Variations in age and metallicity 
are generally correlated, making the dots to align along a narrow 
direction in the plot. The error bars are centered on the mean 
values of the best-fitted Vz and Vt values, and mark the corre- 
sponding 1(T standard deviation. Different labels refer to the CB03 
SSP model (black), and to the M05, CB08, and the t = 1 Gyr 
BC03 models (red). 



the inner and outer radii, Pmin and Pmax, respectively. 
Since our aim is that of characterizing stellar population 
gradients in the ETG population as a whole, we estimate 
the quantities V z and Vt from Eq.[Tl using the mean val- 
ues of the color gradient estimates Vg_x (see Sec. 5). 
We accomplish this in two steps: (i) We parametrize 
the colors available in terms of age and metallicity by 
using different spectrophotometric c odes. We take sim- 

J ile st ellar pop ulations (SSPs) from iBruzual fc Chariot! 
20f)l (BCn3) ■ iMarastonI (|2f)f)5h (M05), and Chariot and 
Bruzual (2008, in preparation; CB08). These SSPs are 
based on different synthesis techniques and have differ- 
ent IMFs. The M05 model uses the fuel consumption ap- 
proach instead of the isochronal synthesis of BC03 and 
CB08. Th e CB08 code implements a new AGB phase 
treatment (|Marigo and Girardil [2007f) . The IMFs are: 
Scab (BC03), Chabrier (M05), and Salpeter (CB08). 
Moreover, we also use a composite stellar population 
model from BC03 having exponential star formation rate 
(SFR) with e- folding time of r = 1 Gyr ®. The models are 
folded with the grizYJHK throughput curves, and flux 
values computed for different values of t and Z. Here, we 

^ We verified that applying the same procedure to the single 
color gradient estimates, and then computing the mean values of 
\/ z ^'iid Vt changes these quantities by less than 0.02dex. 

* Though the determination of absolute values of age and metal- 
licity at a given galaxy radius is beyond the scope of this pa- 
per, we have introduced the model with r = IGyr since it 
is able to match also the ab solute values of central colors in 
ETGs l ILa Barbera et al.|[20M) . 
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use ages spanning from 7.8 to 12.6Gyrs ^, and metallici- 
ties from 0.2 to 2.5^0 for BC03 and CB08 models, and 
from 0.05 to 2Zq for M05. To calculate the color deriva- 
tives with respect to t and Z, for each waveband we fit- 
ted the corresponding flux values with a two-dimensional 
eight order polynomial in log t and log Z. The rms of the 
fits is smaller than O.Olmag for all the bands and models, 
(ii) We solve Eq. [1] in a sense, by minimizing the ex- 
pression 



X' 



.x-XtVt-XzVzY 



(2) 



X 



with respect to Vt and V^. We evaluate the derivatives 
Xf and Xz at a given t and Z values, representing the 
average age and metallicity of the two stellar populations. 

5. RESULTS 

Table 1 lists the statistics oi g — X color gradients. 
We compute the peak value, /i, and the width, a, using 
the bi- weight location estimator (jBeers et al.lll990() tak- 
ing 1000 bootstrap iterations to estimate uncertainties. 
As we can see, the peak value shifts to more negative val- 
ues as we move to longer wavelengths, with a minimum 
fi ~ —0.3 in g — K, and the width becomes progressively 
larger from the r to the K band. Because we have many 
galaxies, the random errors on /i and a are very small. 

Fig-IDshows the mean color gradient, Vg_x, as a func- 
tion of the effective wavelength of the filter X, where the 
trend reflects how stellar population properties vary, on 
average, as a function of radius in ETGs. We adopt the 
procedure described in Sec. [4] to infer the mean age (Vt) 
and metallicity (V^) radial gradients. To this end, we 
minimize Eq. [2] for each of the stellar population models 
described in Sec. HI We estimate the quantities Xt and 
Xz in Eq. dfor t = 10 Gyr i° and Z = Zq. Then, in- 
serting the best-fit values of Vt and Vz into Eq. [Tl we 
derived the color gradient values that best-fit the obser- 
vations. We performed N = 1000 iterations, shifting, 
each time, the observed mean color gradients according 
to the corresponding uncertainties. Each set of best- 
fitted color gradients defines a polygonal line in Fig. [U 
obtained by connecting the best- fitting values of Vg_x- 
The grey region in the figure shows the area occupied 
by the N — 1000 bootstrap solutions. The inset shows 
the values of Vt and V z obtained from the different it- 
erations. Computing the average and the standard de- 
viation of these values, we obtain Vt = 0.13 ± 0.02 and 
Vz = —0.38 ±0.02, respectively. Using a different radial 
range for computing color gradients, with pmin = O.OSre 
to pmax — 2re, docs uot chaugc at all the values of the 
metallicity and age gradients. We also estimated Vt and 
V^ by computing Xt and Xz for different values of t 
and Z, varying t between 9 and 12 Gyrs, and Z between 
3/4 and 2Zq. The inferred absolute value of V^ ranges 
from ^ 0.3 to ~ 0.4 dex, while the value of Vt is always 
positive, varying between ~ 0.01 and ~ 0.25dex 

^ corresponding to a formation redshift of 2 ~ 1.2 and to the age 
of the Universe in the adopted cosmology 

In the adopted cosmology; this age value corresponds to 2 

2.3 

We obtain the minimum value of 0.01 only for the CB08 
model, and only for the lowest metallicity case. For the BC03 
and M05 models, we find Vt>0.1 dex regardless of the t and Z 
values. 



6. DISCUSSION 

We studied mean color gradients in ETGs using an un- 
precedentedly large and homogeneous sample with data 
from optical to NIR wavebands. We find that, on aver- 
age, the main driver of color gradients is a radial varia- 
tion of metallicity, and that a small, positive, age gradi- 
ent seems to be present, on average, in ETGs, imply- 
ing a mildly younger stellar population to the galaxy 
center. The mean metallicity gradient varies between 
Vz = -0.32 ± 0.02 and Vz = -0.425 ± 0.02, depend- 
ing on the stellar-population model adopted to fit the 
color gradients. Previous studies based on much smaller 
samples, found Vz mostly between -0.2 and — 0.3. with 
a typical uncertai n ty of ~ .1 (e.g. iPeletier et al.lll990l: 
Saglia et all l2000l: lldiart et al. 200^ : iLa Barbera et all 
20031 iTa-mura and Ohtan 2003t ide Propris et all 12004 
Wu et al.ll2005ir consistent with what we find here. 
Past literature has shown that age gradients do not 
explain color gradients, although a small positive age 
gradient of Vt ~ 0.1 is still consistent with observa - 
tions (see e.g. lSadia era Dl2000':'La Barb era et~a^l2003^ . 
For instance, iTamura and Ohta (2004) find a value of 
Vt = 0.1 ± 0.14, while iWu et all (|20p), analyzing 
optical-NIR color gradients for 36 nearby ETGs, report 
V = 0.02 ± 0.04. Here, we are able to detect a small 
but significantly positive age gradient in ETGs. The 
presence of younger stars in the center of ETGs is ex- 
pected in most hierarchical formation scenarios. During 
the merger of gas-rich systems, gas dissipates its kinetic 
energy, falling into the galaxy center and forming stars. 
The detection of a positive age gradient also adds 
new insight to the fact that ETGs have smaller ef- 
fective radii at higher redshift when compared t o 
~ (e. ^. J >addi et al. 2005; Trujilb et al,l [200 



Long hetti et alfeoOTl: ICimattr et al.ll2 008: Rcttura et 
i 2009i: Ivan der Wei et al.ll2008t Ivan Dokkum et al.. .20081: 



iBuitrago et all 120081 : ISaracco et al.ll2009f) . Sincethe lu- 
minosity evolution of a stellar population is approxi- 
mately independent of its metallicity, the internal metal- 
licity gradient of ETGs is assumed not to modify the ef- 
fective radius with redshift. However, a difference in the 
formation epoch of two stellar populations corresponds 
to larger and larger differences in luminosity as we ap- 
proach the formation epochs. Our result implies that 
a younger stellar population in the galaxy center would 
brighten more rapidly with redshift than the outer stel- 
lar population, making the profile more concentrated at 
high z, and possibly inverting the observed colour gra- 
dients. To quantify these effects, we take an r^/^ law in 
the same radial range of Pmin to Pmax where we derived 
the internal color gradients (Sec. [3]). At the inner radius 
Pmirn wc assume a stellar population with metallicity 
Zmin and formation redshift z/. Then, for a given ra- 
dius, the profile is evolved with redshift according to the 
luminosity evolution of a SSP whose age and metallicity 
are computed from Zf, Zmin, Vz and Vt. At a given 
redshift, z, the profile is fitted with an r^/* law, yielding 
an effective radius i?e,z- We used three models with dif- 
ferent values of Zf, Vt, and Vz. We fixed Zmin — Zq, 
since varying it in the range of 1/2Zq to 2Zq makes 
i?e,z change only a few percent. All the evolved pro- 
files were well fit by the r^/* law, with an rms smaller 
than ^ Q .Qhmag / arcse(? . Fig. [2] compares the evolu- 
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Fig. 2. — Evolution of the effective radius of ETGs as ex- 
pected from different radial gradients of age and metallicity. Line 
types correspond to different age and metallicity gradient models, 
as shown in the lower— left of the figure (see the text). Circles with 
1(T erro r bars show the ob served ev olution in size fro m previous 
studies llvan der Wei et al.l f2008i) . iSaracco et al.l (|2009|) V 

tion in size, Re,z/ Re,z=Oi from different sources in the 
literature with that expected from the three previously 
mentioned models. Since the data sample different rest- 
frames, from UV to optical wavebands, we constructed 
the models in both the u and g bands. Model A takes 
Zf — 2.8 and Vt — 0.08, an age gradient consistent with 
our findings and such that the stellar population at Pmax 
has an age limited by that of the Universe. Model B 
uses the same as model A, but with Vf = and 
V z = —0.3. Model C has the same metallicity and age 
gradients as Model A for BC03 SSPs, but a younger in- 
ner stellar population with = 1.8. As expected, a 
pure metallicity model (Model B) does not predict any 
significant size evolution with redshift. The evolution is 



mildly stronger in the restframe UV, with Re,z/Re,z=o 
smaller in the u than i n g band. At z ^ 0, the value 
of V„_g is about -0.2 ( IWu et al."2005'). Following the 
method of ISparks and J orgcnscn (1993), this gradient 
implies that effective radii of ETGs increase hy 6 ^ 20% 
from g to u band at 2 0. The waveband dependence 
of Re,z I Re,z=o would tend to reduce and possibly invert 
this trend, producing Vu-g > at high redshift. The 
inversion takes place at a redshift, 27, where the ratio of 
Re,z/ Re,z=o between the g and u bands equals the value 
of S. For both models A and C, zj is very close to 2/, be- 
ing 1.6 and 2.4, respectively. Interestingly, positive color 
gradients in h i gh-red shift ETGs have been detected by 
iFerreras et al.l (|2005D . who found that about one-third 
of field ETGs at 2 ~ 0.7 have blue cores, in contrast 
to only 10% at lower redshift (s ee also iMenanteau et al.l 
[2OOI 12004 [Ferreras et al.l[2009h . 

Fig. [2] shows that a single formation epoch for the inner 
stellar population of ETGs does not reproduce the ob- 
servations. Model A predicts a size evolution of ~ 80% 
at z ~ 2.4, in good agreement with the data at z > 2. 
However, the evolution is too shallow for lower redshifts, 
where a smaller z/ is required to fit the data (model C). 
This result could be explained in a picture whereby ETGs 
form at different redshifts by gas-rich mergers, and the 
redshift-size trend is due to the increa se in the amount of 
dissip ation in mergers with redshift (jKhochfar and Sil3 
I2OO6I ). Fig. [2] suggests that the central bursts of star for- 
mation take place at different redshifts, above z 1.5. 
However, the present analysis does not exclude a picture 
where the younger stars in the center form at very high 
redshift in most galaxies, with the color gradients ac- 
counting for only some part of the redshift-size relation. 

We thank G. Djorgovski, I. Ferreras, R. Gal, and 
P. Sarracco for several comments that helped to im- 
prove this paper. We thank the referee for help- 
ful su ggestions. We used data fro m the SDSS 
( http: / / www.sdss.org/ collab oration/ credits .html]) . This 
work is based on data obtained as part of the UKIRT 
Infrared Deep Sky Survey. 
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